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ABSTRACT 

Expressions for the excitation cross section of the highly excited states of 
the hydrogenlike atoms by fast charged particles have been derived in the dipole 
approximation of the semi classical impact parameter and the Bom approxima- 
tions, making use of a formula for the asymptotic expansion of the oscillator 
strength of the hydrogenlike atoms given by Menzel. When only the leading term 
in the asymptotic expansion is retained, the expression for the cross section 
becomes identical to the expression obtained by the method of the classical col- 
lision and correspondence principle given by Percival and Richards. Compari- 
sons are made between the Bethe coefficients obtained here and the Bethe co- 
efficients of the Bora approximation for transitions where the Bom calculation 
is available. Satisfactory agreement is obtained only for n -* n + 1 transitions, 
with n the principal quantum number of the excited state. 
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I. INTRODUCTION 


A knowledge of accurate values for the excitation cross section of highly 
excited states of the hydrogenlike atoms by electron impact is necessary in 
order to understand the conditions under which the observed radio frequency 
transitions between the highly excited states of these atoms in the interstellar 
ionized region, or similar processes in the solar corona, take places. Using 
the classical collisions and the correspondence principle Percival and Richards 
(1970) give an expression for the excitation cross section valid for high impact 
energies. Similar expression is given by Presnyakov and Urnov (1970) making 
use of the impact parameter approximation and an atomic model based on an 
harmonic oscillator. 

The work that is being reported here gives solution of the same problem in 
the dipole approximation of the semi classical impact parameter and Bom ap- 
proximations, making use of the asionptotic expansion of the hydrogenlike oscil- 
lator strength given by Menzel (1968, 1969). The results of the calculation are 
checked against the results of the Born approximation. 

In collisional excitation the Bom approximation becomes valid when the 
incident energy is some 50 times larger than the excitation energy. The Born 
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cross section then should not differ from the actual cross section by more than 


a few percent. In low impact energy excitation of highly excited states of atoms 
another consideration should be taken into account. Due to the closeness of 
their energy levels, many excitation channels are strongly coupled to each other, 
and validity of a two channel approximation such as the Born approximation be- 
comes questionable (Seaton, private communication). 

Nevertheless, here we assume the validity of the Born approximation at high 
impact energies, and compare the results of the dipole approximation, itself an 
approximation to the Born approximation, with the results of the full Born ap- 
proximation. Instead of comparing the cross sectional values for the two approx- 
imations, it is more convenient and more useful to compare the corresponding 
Bethe coefficients which are independent of the mass, charge, and energy of the 
projectile. In Table I this comparison has been made. 

We have also shown that when the leading term in the expansion of the oscil- 
lator strength is retained, the result of the impact parameter dipole approxima- 
tion becomes identical to the result of the classical collision and correspondence 
principle method given by Percival and Richards. 

II. FORMULATION AND RESULTS 

Let n be the initial and n' the final principal quantum numbers of a hydrogen- 
like atom with atomic number Z in collision with a charged particle with charge 
Z'e, e being the absolute value of the electronic charge. According to the dipole 
form of the semi classical impact parameter method the excitation cross section 
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between these two levels is given by (Seaton, 1962) 


p(n - n' ) bdb (1) 

HI i n 

where b is the impact parameter and 


o-(n -* n' ) = 2 tt 
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(2) 


with r the position vector of the bound electron, M, E, and V the reduced mass, 
the relative energy, and the relative velocity of the system, m^ the electronic 
mass, and AE the excitation energy. The function f(x) is defined by 

f(x) =x2[k,2(^) + K^^(x)] (3) 


where Kq(x) and K^(x) are modified Bessel functions. 

The matrix in (2) represents the expectation value of r between sub-levels 
of n and n', averaged with respect to the initial states and summed with respect 
to the final states, namely 
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where 'i'i' and mm' are the azimuthal and magnetic quantum numbers. 


4 



Introducing s by s = n' - n then for both n and n' large and s small the 


matrix in (4) is given by (Menzel 1968, 1969) 


l^n' I ^|n>l ^ = 




X _L J ^ (s) J (s) g(s), g(s) = 1 

3s^ n2 


(5) 


where a = n^ atomic radius before excitation and a.^ the Bohr radius. 

t 

J^(s) is the Bessel function of equal order and argument and (s) is the first 
derivative with respect to the argument. A(s)are constants of order of unity given 
by Menzel (1969). Substitution of (5) into (2) gives 


V Ik *7 t 2 ^ 

p(n-n)=Z ^ — 
tn . 




To evaluate the cross section given by (1) we choose b^in to be equal to a. 

A different choice for b . and of the same order of magnitude as a makes little 

min 

difference in the high energy values of the cross section. Making use of a formula 
for an integral over the squared of the Bessel’s functions (Morse and Feshbach, 
1953) we obtain the following expression for the cross section 


(n — n')_ 16Z' MRy ^4 A 


3Z2m E 

e 


j.(s)j:(s) 



By putting g(s) - 1 Equation (7) becomes identical to the expression for the 
cross section given by Percival and Richards (the right hand side of their ex- 
pression should be multiplied by a factor of n'*). Since the impact parameter 
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method is valid when the incident momentum is much larger than the momentum 
of the orbiting electrons, and the dipole approximation is valid when the momen- 
tum transferred to the atom is small, we conclude that the results of Percival 
and Richards are also valid under the same conditions and provided n> > 1. 

Making use of the relation 


^aAE^ _ M AE^ 

/ 4:n^a^RyE 


(8) 


and the validity criteria for the impact parameter method and the dipole approxi- 
mation it can be seen that aAE/(hV) <d always, except for the rare occasion 
that the incident velocity is much less than the orbital velocity. Then an expan- 
sion of the modified Bessel functions in (7) about the origin, as has been done by 
Percival and Richards, is permissible and will yield 


cr(n“*n')_ 8 Z'^MRy ^4 


17 a" 


3 Z^m E 


1 +^ 
n 


J,(S) 13(5) 


g(s)'tn [5.04 


m^a|R^' 
M a 2 AE 2 


•(9) 


It should be emphasized that for the range of the incident energy that (7) is valid, 
(9) can be used instead (Percival and Richards have a factor of 4.5 instead of the 
factor of 5.04 in (9). The reason for this discrepancy is not imderstood). 

To show more clearly the dependence of the cross section on the impact 
energy it is convenient to write (9) in the form 
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(10a) 


^(n - n') _ Z'^ M Ry 


77 a' 


Z2m E 


m E 

A(n,n')'tn( — ) +B(n,n') 

Z2 Ry> 


where A(n, n') and B(n, n ') are pure numbers depending on n and n’ only and 
are given by 


A(n.n') = I (l s"" J,(s) J'(s) g(s), 

3 V n ' 


B(n , n ' ) = A(n , n ' ) -tn [5. 04(l-n^/n'^) ^]. (10b) 


A numerical table for M(s) = 4/3 s~^ Js(®) given by Menzel (1969). 

Similarly, the numerical values of g(s) defined by (5) can be obtained using this 
reference. 

To drive the dipole form of the Bom approximation we start with the ex- 
pression for the Bora approximation given by 


cr®(n n') 


8tt 

E/Ry 



1 '(n' 



( 11 ) 


with kj the wave number of the relative motion before collision related to E by 


(M/m^) (E/Ry), 


( 12 ) 


and the wave number of the relative motion after collision related to k^ by 


k ^ (M/m^) (AE/Ry) , (13) 
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andirq the momentum transfer between the particles. The matrix in (11) in 
analogy to (4) is defined by 

't'm' 

with'tm^'m' defined before. 

When a^^ >> (M/m^ ) (AE/Ry) we have from (12), up to the first order 

terms, 

kj - = [M/(2m^ ag^kj)] AE/Ry, kj+k2=2k^. (15) 

The radial part of the integral <^n' |exp (iq • r) | n)> decreases exponentially 

when 1 r | becomes greater than the atomic radius a = n^ Following Bethe 

(1930) we expand exp (iq • r) under the integral sign and introduce a parameter 

q such that 
0 



E Kn'^' 


iq 
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( 14 ) 


[M/(2m^aQ2k ^)] AE/Ry « « Z/(n^ap), (16) 


then (11) can be written in the form 
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q^ 


( 17 ) 
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In deriving (17) the range of integration in (11) has been divided into two 
parts; the first part running from kj - k 2 to and the second from to k^ + 
k^. Due to (16), only the dipole term in the expansion of exp (i^ • r) in the first 
integral has non-vanishing values. Integrating this term and using (15), we find 
the first term in the bracket of (17). The second integral cannot be integrated 
easily, but in any case the cross section is independent of q^ . From inequality 
(16) it is evident that at high impact energies we can choose q ^ independent of 
the impact energy. With the help of (5) and (12) Equation (17) can then be 
written in the following form 


cr^(n - n') _ Ry 


77 a; 


Z2 m E 


A®(n , n ') 'in 


'm E 
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^MZ^Ryy 


+ B®(n, n') 


AB(n,n')^^n^ (l + - s"' J /s) J ^(s) g(s), 


(18a) 


B®(n , n ' ) - A®(n, n' ) 'in 


4 q2 


(1 - nVn'") 


2n2 


2k, 


Kn' 


iq 


ln>12 Sa, 


n >> 1 (18b) 

Equations (17) and (18) have been derived within the validity of the Born 
approximation. It is seen tliat A® (n, n ') is identical to A(n, n') given by (10b). 
If we now make the approximation that q^ ^ l/a instead of the second inequality 
in (16), that is q^ is inversely proportional to the atomic radius, then (17) will 
be an approximation to the Bom approximation in the sense that in its first 
term on the right hand side we have neglected contribution of all the multipoles 
higher than the dipole. From Equation (18) it can be seen that the choice 
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= 1/a has no effect on the value of A® (n, n'), but it lowers the value of 
(n, n') from its actual value. This is due to the neglect of the higher multi- 
poles in the first term in the expression for (n, n*). In Table I where we 
compare the Bethe coefficients derived here with the Bethe coefficients of the 
Born approximation we find this in fact to be the case. 

For q^ = 1/a the value of the integral in the expression for B® (n, n') is 
small. This is due to the fact that when q becomes several times larger than 
the atomic radius, because of the factor exp (i q - r) the integrand of 
</ n' 1 exp (i q • r) 1 n^ will oscillate rapidly in its range, assuming positive 
and negative values, and the integral will be small. With qQ= 1/a and the neglect 
of the integral, B® (n, n') becomes the same as B(n, n ') given by (10b) except 
for an additional factor of 1.26 in the argument of the logarithm in (10b). 

The Bethe coefficients in the Bom approximation. A® (n, n') and B® (n, n'), 
for n = 1 - 9 and many values of n' in the range 2-20 are given by Omidvar 
(1969). In Table I the ratios A(n, n')/A® (n, n') and B(n, n')/B® (n, n'), where 
A(n, n ') and B(n, n ') are given by (10b), are given for a range of n from 5 to 9 
and a range of s from 1 to 4. As is seen the A ratios range from 0.81 to 0.95. 
The reason that these ratios are different from unity is the use of the asymptotic 
expansion of the oscillator strength for calculation of A(n, n') instead of the 
oscillator strength itself. The B ratios are, however, quite different from 
unity except for the n n + 1 transitions. The reason for the smallness of these 
ratios is the approximation q^ = 1/a. Physically this approximation is equivalent 
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to an arbitrary cut off of the large momentum transfers to the atom, and a 
resulting smaller cross section. 

We conclude that the dipole approximation of the semi classical impact 

r 

parameter or Born approximation gives as expected the A ® (n, n ) coefficient 
correctly, but not the B®(n, n') coefficient. Contribution to the cross section 
coming from the B®(n, n') coefficient even at moderately high impact energies 
is appreciable. Then an adequate description of the cross section by the dipole 
approximation for moderately high impact energies, with a possible exception 
of n- n + 1 transitions, can not be given. It has been shown that the method of 
classical collision with correspondence principle is equivalent to the dipole 
approximation. 
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Table I 

Ratios of the Bethe Coefficients. For each n and s = n' - n two numbers are 
specified. The top number is A (n, n')/A® (n, n ') and the bottom number is 
B(n, n')/B^ (n, n'). A(n, n') and B(n, n') are the Bethe coefficients in the semi 
classical impact parameter, and A® (n, n') and B® (n, n') are the corresponding 


coefficients in the Born approximation. 


n 

s\. 

5 

6 

7 

8 

9 

+1 

0.9176 

0.9292 

0.9378 

0.9451 

0.9510 

0.8015 

0.8015 

0.8033 

0.8069 

0.8132 

+2 

0.8646 

0.8800 

0.8922 

0.9026 


0.4064 

0.4111 

0.4174 

0.4250 


+3 

0.8309 

0.8460 

0.8594 

0.8713 


0.2741 

0.2739 

0.2768 

0.2913 


+4 

0.8064 

0.8200 

0.8326 



0.2117 

0.2081 

0.2079 















































